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Staggered fermions and their 0(a) improvements 
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a Columbia University, Department of Physics, New York, NY 10027 

Expanding upon the arguments of Sharpe, we explicitly implement the Symanzik improvement program demon- 
strating the absence of order a terms in the staggered fermion action. We propose a general program to improve 
fermion operators to remove 0(a) corrections from their matrix elements, and demonstrate this program for the 
examples of matrix elements of fermion bilinears and Bk- We also determine the additional operators which must 
be added to improve the standard staggered fermion currents. 



1. INTRODUCTION 

In order to reduce the systematic errors com- 
ing from finite lattice spacing, we must improve 
both the action and the lattice operators. Accord- 
ing to the improvement program of Symanzik, we 
find all dimension-5 operators which are invariant 
under the lattice symmetry group and add them 
to the original action to reduce all 0(a) correc- 
tions. If there exists no such dimcnsion-5 opera- 
tor, then the action is already good to 0(a 2 ). On 
the other hand, even though the action is accurate 
to 0(a 2 ), the matrix elements generally may still 
have 0(a) errors, and hence, we have to improve 
the operators themselves. 

In this paper, I will expand upon the argu- 
ments given by Sharpe M to prove that there are 
no O(a) terms which can be added to the stag- 
gered fermion action, and hence the action is al- 
ready accurate to 0(a 2 ). Then I will propose a set 
of improved fermion field variables and use them 
to construct the fermion operators to reduce the 
0(a) corrections to their matrix elements. We ap- 
ply this program to the case of (0|s^54C?| AT ) and 
Bk as examples. We find the former does have 
0(a) corrections while the latter does not. We 
will also determine the additional operators that 
must be added to improve the standard staggered 
fermion currents to define operators whose matrix 
elements are accurate to 0(a 2 ). 



2. STAGGERED FERMION ACTION IS 
ACCURATE TO 0(a 2 ) 

What we will prove here is that there exists 
no dimension-5 operator which is invariant under 
the lattice symmetry group (rotations, axis rever- 
sal, translation, U(l) ®U(\), charge conjugation, 
etc), and therefore, the staggered fermion action 
is already good to 0(a 2 ). 

We rewrite the staggered fermion action as 

S/=xE(t^®W* + "»]x, (i) 

and consider all dimension-5 operators which 
have the general form Xls ® £,fJ(D)x where / 
is a homogenous real polynomial of degree 2. 

Invariance under E/a(1) requires that S + F is 
odd, so only the following combinations of S ® F 
are valid: 



(-^,75,7^.75^) ® (£a,6sa), 

(7m: 75 M ) ® (I,£,5,£,\t,^5\t)- 



(2) 
(3) 



The axis reversal invariance will further limit 
the number of terms to be four: 



75 ® 

T&[7m>7*]®&(£m+&0[ £, /*>A/], 

75 bv» iA ® - 60Pm. d »}- 



(4) 
(5) 
(6) 
(7) 



The rotational invariance will further eliminate 
the term in Eq. (Q) but allows the remaining three 
terms Eq. (| - g) . 
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Finally, none of the terms listed in Eq.(||- 0) are 
invariant under lattice translation! So, we con- 
clude that there is no dimension-5 fermion opera- 
tor which is invariant under the lattice symmetry 
group, and therefore no dimension-5 operator can 
be added to the staggered fermion action. 

3. IMPROVE STAGGERED FERMION 
FIELDS 



Following ref . p| , we define the hypercubic fields 



as 



1 A 



(8) 



(9) 



Then the classical continuum limit of the stag- 
gered fermion action can be written as: 

S F [q,q}^ I ^«(2/){X(7„<»-0ab£>m 

J y ab 



iSab - a[y,(75 ® £,5h)abD 



' 4 



®&(£n+&))ABF^(y)}}q(y) 
+0(a 2 ). (10) 

where = d^ + igA^ is the continuum covariant 
derivative. At first sight, the action contains or- 
der a terms. Likewise, it is clear that the fermion 
propagator for the hypercubic fields q and q devi- 
ates from the continuum propagator by terms of 
order a. However, if we introduce the following 
improved field variables 



x I A (y) = (i-aJ2^D^)xA(y), 
x A {v) = XA{y){l - aJ2 A » S), 



(11) 



(12) 



and replace Xi X m Eq.(|8| ||) by x 1 an d X 1 > we 
will reduce the finite lattice spacing corrections 
from O(a) to 0{a 2 ). 



Using the new fermion fields, we can construct 
improved fermion operators. For example, the 
improved fermion bilinears have the following 
form: 



Xa{v){1s ®£,F) AB XB{y) 



XA{y){is ® £f) ABXB(y) 



AB 



■{lhuS®£,5uF) AB )XB{y)] 



'^2xA{y)[{lhuS ® &i/f) 



AB 



-{lS5u ® &5v) AB ]D v XB{y) 



+0(a 2 ). 



4. APPLICATIONS 



(13) 



Here, we apply this improvement program to 
the calculation of matrix element (0\sj^4d\K °), 
the calculation of Bx , and the improvement of 
lattice currents. 

4.1. (0|s 754 d|^ ) 

The axial current used in the (Landau gauge) 
numerical simulation is: 



A »{y) = * A (75m® &) A b Xb(v)- 



(14) 



AB 



From the continuum expression 
P(t) cont = (0\Ai{t) cont \K Q ) 



= V2f. 



K m k e 



-m k \t\ 



we define, on the lattice, 

p(t) = {o\yM(x,t)\K°), 



(15) 



(16) 



and put the wall source that creates the K° 
on the time slice at t = 0. If we don't 
consider 0(g 2 a) terms, we can take only the 
term -f 9^( 7 54 ® & ) AB Xb] in Eq.© 
because other terms contribute zero "flavor" trace 
at the tree level. So, we have 



P(t) Im 'P = P(t) - -d^P(t), 



(17) 
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and this simple example show that, if we do not 
consider O(gQa) corrections, the improved opera- 
tor gives the well-known t + | interpretation: 



At*< I (y) = Af(y) 



-m k \t+l/2\ 



P{t) = V2f K m k e 

4.2. B K 

The formula for calculating Bk is: 

M K 



B K = 



3 JVl K 



(18) 



(19) 



The improved numerator is (omitting the 0(g^a) 
terms): 



M 



Imp 
K 



M 



K 



-diM 



K 



0(a 2 



(20) 



Since J\4x(t) is computed from a plateau (i.e. 
time independent) within the statistical error, 
there is no O(a) corrections to the numerator. If 
we calculate the quantity (0|s7475G?|i4T°) in both 
forward and backward directions and multiply 
them to get the denominator, then there will have 
no O(a) correction even if no attention is paid to 
an accurate definition of fx. Hence, we showed 
that there is neither 0(a) nor 0(ag 2n log" a) cor- 
rections to Bk- Sharpe ffl has examined this 
question in greater detail and argued that in fact 
there are no corrections of 0(g 2 l n a) also. How- 
ever, if we calculated the denominator only in one 
time direction, there would be an error of order 
of 0(mka). 



4.3. Renormalization of lattice currents 

The lattice currents can be written as||: 



Jfatt = XI J ®£fX, 



(21) 



Using the method developed in this paper, we can 
explicitly determine the improved currents accu- 
rate to 0(a 2 ). For example, the conserved vector 
current and axial vector current corresponding to 
the Uv(l) £S> Ua{1) lattice symmetry can be writ- 
ten as follows: 

Vfa) = V,(y) - 9u[xa(^I)abXb] 



~2 9 m[xa(7m® I)abXb] 
^2dv[XA(lZ,[^v] ®&v)abXb] 



"o X! 9 "Ixa(7h5 ® £,s)abXb] 



- 2^[xa(7^5 <B> ^)abXb] 



+ ^ d ^XA(l[^] ® £u)abXb] 
+0(a 2 ). (23) 



+0(a 2 ), 



(22) 



The effect of the second term on the right hand 
side is to shift the position y, labeling the cur- 
rent, from the corner to the center of the hyper- 
cube. The third term whose effect is to shift in 
the jti's direction occurs here because the currents 
are non-local operators which involve an overlap 
between two nearest hypercubes. The forth term 
is a mixing of a different spin-flavor operator and 
is necessary to remove all order a effects from a 
general matrix element. 

5. CONCLUSIONS 

The staggered fermion action is accurate to 
0(a 2 ); In order to reduce 0(a) corrections in 
a matrix element, we must use the improved 
fermion fields which we have proposed; Our im- 
proved program derived the well-known t + ^ in- 
terpretation for the one hypercubic fermion bi- 
linears at the tree-level; We demonstrated Bpc is 
good to 0(g 2 a) at the tree-level; We proposed 
terms which must to included in the lattice cur- 
rents in order to reduce 0{a) corrections. 
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